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A new class of static plane symmetric solution of Einstein field equation, which is judged as 
the source of Taub solution, was presented in our previous work. In this letter the properties of 
geodesies of this solution are explored. It is found that this solution can be an appropriate simulation 
to the field of a uniformly accelerated observer in Newton mechanics. The essence of the source is 
investigated. A phantom with dust and photon is suggested as the substance of the source matter. 
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I. INTRODUCTION 

We successfully found a new class of static plane symmetric solution of Einstein field equation sourced by a perfect 
fluid in [1]. A very general 4-parameter solution was presented and then a reduced 2-parameter family was studied in 
detail 

This solution puzzles out a long standing problem: the source of Taub solution 2]. A no-go result is shown in 
which says that a perfect fluid cannot bound a vacuum in a space with plane symmetry unless the boundary 
condition of the continuity of the derivatives of the metric tensor is violated. This means that there does not exist 
a matter source which can perfectly match to vacuum Taub solution. A requirement is imposed in the proof of 
this no-go theorem, that is, the pressure is positive throughout the slab for the solution. However, it is shown that 
generally speaking any singularity free source with reflective symmetry for plane symmetric vacuum space must violate 
' dominant energy condition (DEC)[3. In fact, we find that DEC is always violated in the solution in no matter 
^ , what values the parameters are taken. Further, we have found that there is a configuration in the class of [l| which 
can perfectly match to vacuum Taub solution, and it is naturally to be identified as the source of Taub solution. 
To reveal more physical significance of the solution, we need to study the motions of massive particles (in the 
i following context called particles) and massless particles (in the following context called photons). We will study the 
• ■ geodesies of particles and photons respectively. Although the energy-momentum form of the solution is obtained, we 
i; have no idea what it really is. We will study the possibilities for various candidates, including dust, photon, scalar, 
Qv^ ' and phantom. 

, This letter is organized as follows: In the next section we present a brief review of the 2-parameter family in [1] as 
IL* ' our starting point for further discussions. In section III we investigate the geodesies in this space. In section IV we 
. ^ ' study the substance of the source. A summary and some discussions are presented in section VI. 

P_ , II. THE SOLUTION 

The plane symmetric solution sourced by perfect fluid, which is presented in fl[, reads, 

ds^ = -e^-^dt^ + dz^ -f e^^'''+'""\dx^ + dy^). (1) 

There are two free parameters a, b in the above solution. With the above metric, p{z) and p{z) are respectively given 
by 

p{z) = -a2(3 + 86e°" + 352e2"^), (2) 
p{z) = a2(3 + 45e°^ + 52g2az-)^ 
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where we have set 8ttG = 1, and G, as usual, denotes the Newton gravitational constant. ([2|) and ([3]) determine an 
implicit function p = p{p). It can be treated as the most general form of the source of plane symmetric space in (5|]. 
If we require it to match to a vacuum space, p should vanish naturally at some distance zq from the ground z = 0, 

pU=zo = 0. (4) 

Under the conditions 

b = -3e-°^o, (5) 

azQ = -1/3, (6) 

the above solution can perfectly match to Taub solution (Here, we present a minor correction to [l|, which does not 
change any conclusion of [l|). As we have shown in [l|, the metric (P) degenerates to AdS (anti-de Sitter), 

ds^ =dz^ + e^''''{-df + dx^ +dy^), (7) 

when b goes to zero. 

III. GEODESICS 

Free falling is natural state of a particle. Mathematically, its orbit is time-like or null geodesies. In fact, any 
test particle with the same initial position and velocity goes along the same orbit: It is just the (weak) equivalence 
principle. Fundamentally, we can derive the geodesies by the self-parallel equations, 

d^x'' „ da;" dxf^ 

where x'^ = {t,z,x,y), T denotes the afhne connection of metric ([1]), r is an afhne parameter along the geodesic. 
We can calculate the afhne connections routinely and then solve the four coupled second order equations. Generally 
speaking it is not an easy work. 

Alternatively, we advance a different way by applying the properties of Killing vectors. The velocity normalization 
imposes the constraint along a curve, 

,dd. f dt\^ / dz\'^ / dx\^ f dy\^ 

'^d-r^ d-r^ = ^" U j + U j + U ) + [d-r ) = 

where g is the metric tensor whose components are given by ([l}, rj = —1,0,1 for time- like, null, space- like curve, 
respectively. For a time-like curve, r denotes the proper time; For a null curve, r represents an affine parameter; and 
for a space-like curve, r stands for the length parameter. The inner product of a Killing vector ^ and the tangent 
vector r = ^ of a geodesies is a constant along the geodesic. This result can be proved as follows, 

VT(g(r, 0) = 5(VtT, + .9(T, VtO = 0, (10) 

where V is the derivative operator consisting with the metric g. The three Killing vectors ■§[, and ^ yield 

.(T,|)^..|^-i., (11) 
9{T,^)^g..^ = Pu (12) 

respectively. Here, E, Pi,P2 are three constants. Physically, for a time- like geodesic E denotes the energy of a unit 
mass particle moving along it, Pi,P2 are the momentum of unit mass particle along x and y direction, respectively. 
From now on we only consider time-like and null curves since the particles of realistic matter can run only along such 
curves. 
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First we consider the motion along z direction, that is x =constant, y =constant. In this case we should solve the 
equation set (O and pT|) . Note that the geodesies along z direction will be the same as ones of an AdS, since the 
equations to determine it are exactly the same as that of AdS, say, Q and (fTTjl . The reason is that the AdS metric 
([7]) permits the Killing vector ^ and the tangent of the geodesic also needs to obey the normalization equation ([5]), 
where gu and g^z are the same with ([T]). For time-like geodesic, the result reads 

z = - [In £" + In I sin(aT + ci) |] , (14) 



t = cot(aT + ci) + C2, (15) 

where ci, C2 are integration constants. It is clear there is a maximum value of z. For such an observer, the bigger 
value of i?, the farther it can go. There is an apparent problem for this observer, z and t will be divergent when 
AT + ci — > riTT, where n is an integer. This indicates that z,t are not suitable coordinates at the neighborhood of 
or + ci — nix. The case is very similar to the case of Schwarzschild coordinate t. The Schwarzschild coordinate t of a 
free falling observer will be divergent when it goes to the event horizon. To find a well-posed coordinate is not very 
difRcult. For example, we set a new coordinate / to replace 

1= i + (^'Ic^)2 =«gsin(2(T-fci)), (16) 



and a new coordinate J to replace z. 



J = ae^ = £^1 sin(ar -I- ci)|. (17) 



It is easy to see that for any r, both / and J keep finite. 

Now we consider its 3- velocity measured by t. The 3-velocity is defined as 



where i runs from 1 to 3. Here only the velocity in z direction does not vanish, whose magnitude reads 

dz ^ dt 
dT dr 



= [9Uv,vr/'^^A-\ (19) 



By using (fT4)) and p5|) . we arrive at, 



Expanding v about t ~ C2, we obtain. 



1 + (C2 - t)^a^E^ ■ 



(20) 



a^E^it - C2) + a*E\t - c^f + O {{h ~ cf) . (21) 



In the region around t — C2, constant acceleration is a fairly good approximation since there is no (t — 02)^ term 
in the expansion series. In (j21[) . v is independent of z,x,y. In such a sense, the global effect of gravitation, not 
only the local effect, can be simulated by a field of a uniformly accelerating observer. One may doubt that this is 
a physical conclusion since a coordinate system can be rather arbitrary, that is, if we change a coordinate system 
whether this conclusion remains valid. We present a brief analysis of this point. ^ is the time-like Killing vector 
which denotes that the metric is static. Hence the coordinate i, as its integral curve, is unique (up to a constant 
factor). Similarly, the coordinates x^y are unique (up to constant factors). The coordinate z is orthogonal to all of 
them in the 4-dimensional space-time with g^z = 1. Hence it is unique too. So our conclusion is physical. A research 
on plane symmetric solutions in order to find the best simulation to general relativity of the Newtonian infinite plane 
was presented in Q. For some interesting researches on the equivalence between a uniformly accelerating reference 
frame and the gravitational field, see Q. 
For null curves, ?7 = 0, we reach, 

z^-\n{aET + C3.), (22) 
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t=- 2F2\ (23) 

where C3, C4 are integration constants. 

For the most general case in which x and y are not constants, we have to solve the associated equation set (jH [TTl 
[T2I [T3|) . In this case a particle goes along a different geodesic since g^x and (?yj, in ^ become different from that 
of ([7]). We can derive some physical conclusions directly from the equations without analytical solutions. Since the 
space is inhomogeneous in z direction, the geodesic along which z =constant does not exist. So any initial particles 
moving in a; or y direction also must move in z direction. However, a test particle will go through the same proper 
distance in the same proper time along x or y direction in such motions, that is, so it moves uniformly in the x — y 
plane. For example, along x direction 



7 = (24) 

or 

where we have used (fT2|). For the case of null geodesies, we have the same conclusion, but we need to change the 
"proper time" to "afhne parameter" . 

For a time-like geodesic, it is very difhcult to find an analytical solution in the case of x or y is not constant. For 
a null geodesic, we find an analytical solution in integral form for y =constant. In fact it is the most general case 
because the plane symmetry exists on the x — y plane. It is a solution of equation set ^ llll I12p with 77 = and 
y ^constant, which can be written as follows, 

/dz 
, 25 



where /3 — Pi/E is a constant which is smaller than 1, and 

dz 



-26e 



(26) 



0(1-^2)1/2 ^2)3/2 2a(l-/32)5/2 



Here we directly write the null geodesic by using the coordinates rather than r, for now r has no significant physical 
sense. Our solution ([1]) can be treated as a generalization of thus we expect the geodesies will degenerate to that 
of AdS in the limit 6 0. We have shown that the motion equation along z direction is just the same as that of AdS, 
since b does not appear in that equation. Now we expand (|25p and (|26p into series about 6 = 0, respectively. 



and 



where C5 and cg are two integration constants. The leading terms in (j27( I28 p are the corresponding terms in the case 
of AdS. The higher order terms denote the corrections to AdS. 

IV. THE ESSENCE OF THE SOURCE MATTER 

Although we get the exact form of p and p in ([2]) and ([3]), we do not know what it really is. We will investigate 
some different candidates for the source matters. 

First, it can not be pure dust or photon, since p ^ and p ^ ^p. 

Next, we consider the source matters are consisting of dust and photon. In this case, 

P = Pdust + Pphoton, (29) 
P = Pphoton- (30) 



Pphoton and Pphoton are related by 



J'photon — ^/^Photon- {^^) 
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By using ([2]) and ([3]), we derive 



(32) 



and 



Pphoton = 3a' (3 + Abe^' + b'e""') 
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(33) 



Pdust = -^-[2 + —be- + b^e^''^). (34) 
For the special configuration whicli can perfectly matcli to Taub space, b = — 3e^/'^, azp = —1/3, pdust becomes. 



Pdust — —12a 



^ 4 ' ^ 16 



< 0. 



(35) 



Thus, wc find an unusual property of the source matter: pdust is always negative. Though it can be realized by some 
quantum effects Q, in classical level, we may need to seek more natural matter for the source. 

Third, we prove that the source is not a scalar field. In coordinates ([T]) the components of the energy-momentum 
of a scalar read, 



(36) 



(37) 
(38) 



where 1^(0) denotes the potential of the scalar field. Einstein equation requires 



(39) 



(40) 



(gH) yields. 



and 



= e-('^^+'"=")9,0. 



(41) 



(42) 



(43) 



From now on we only consider the static solution. Under this assumption, we derive the general solution of (|42p and 

(1131), 



= C 



-(az+be") _ bEii-be"'') 



where C is an integration constant, and Ei denotes the exponential integral function, which is defined as 

^ f°° e"'' 
Ei{u) = — ds . 



From (|40| and ((4T|) . we obtain another equation that cp must satisfy. 



(44) 



(45) 



(46) 
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However, does not satisfy the above equation p6)) . To show this, we define 



Direct calculation presents, 



(47) 



(48) 



We see that F depends on z and hence cannot be identical to zero. So it is impossible that a scalar plays the role of 
the source matter. 

Observing the energy momentum of metric H]) carefully, we find that the energy density is prone to be a negative 
number. Especially in the case of matching a Taub solution it is negative. It is natural since this solution is a 
generalization of AdS. Thus, we should consider sources dominated by phantom fields. The phantom has been widely 
investigated in cosmology after the discovery of cosmic acceleration. It can violate the domination energy condition, 
which is a requirement of a reasonable source for Taub solution And it just can be realized in metric (|T]) Thus 
it would be proper to investigate the possibility of a phantom source. The discussions are parallel to the case of a 
scalar. The energy-momentum for a phantom ^ can be written as follows. 



(49) 



(50) 



(51) 



y 



-(-Ig'^dti^M - Ig^'dMib - ig^^a.VS.V + Ig'^'dyi^dy^) - (7(V), 



where U{'ijj) denotes the potential of the phantom field. Similarly to the case of a scalar, we derive 



In the context of a phantom, ([46]) becomes 



v2„2„2az 



And defining 



we find 



(52) 



(53) 



(54) 



(55) 



(56) 



where D is an integration constant. Therefore is a function of z and a single phantom can not behave as the source 
either. 

We see that the source matter may be much complicated than we expected. It may be composed by several different 
components. Assuming that every component behaves as perfect fluid, that is, it is isotropic. 



Pix Piy Piz t 



(57) 



where pi denotes the pressure of the i-th component. The constraint equations for ■(/;, similar to that of (j) (j42p and 
(|43)) . are still valid. The solution for general solution for ■0 is the same as above. 



^ = D L-(a^+be") _ ^ ^ bEii-be"') 



(58) 



7 



Now we try to derive the potential U{i/j) through the motion equation of the phantom. Our strategy is as follows: One 
first derives the potential through the motion equation of the phantom. Second one substitutes the potential to the 
energy momentum form of the phantom (|49ll52p to obtain the density of pressure. And then one compares the energy 
momentum form of the metric ([1]) with the energy momentum form of the scalar to find which are the other necessary 
components of the source. 

The motion equation of the scalar reads, 



v/-det(.9) 



(59) 



where det{g) marks the determinate of the components of the metric ([T]). Substituting (|58p into above equation, we 
derive 



U = -a^D^ 



{be"'' - i)e-2('^^+f«="') + 2b^Ei{~2be''' 



where Di is an integration constant. Then, the density and pressure of the scalar read. 



D 



1) 



(60) 



(61) 



J^phaiitom — 2^ ^ 



2^2„-2(oz+2be" 



-1 + 266°^ + 4&2e2("^+2'"=")i;i(-2&e° 



A. 



(62) 



If the source is composed of other components except for the phantom, the total density pot and pressure pot of other 
components can be written as. 



p„t = -a2(3 + 86e''^ + ih'e^''^) + ia2£,2g-2(a2+26e") ^ c^^^az _^ 4^2g2(az+2fce")^-(_2^ga2 

2 



(63) 



= 0^(3 + 46e"- + 62e2-) _ ^a2^2g-2(a.+26e") 



-1 + 26e"^ + 462e2(«^+2''e''")^;j(_26e"^) 



(64) 



We see that the pressure of the components other than the phantom does not vanish. Therefore, a single fluid of 
dust can not play the role. The pressure of the other components is not 1/3 of the density, hence a single fluid of 
photon can not play the role either. We consider a source composed of a phantom, dust, and photon. The pressure 
and density of the photon should be. 



„ „2^Q I AU^az I i,2„2az\ ^ 2 n2 -2(az+26e"^) 

Pphoton = a (3 + 40e + o e ) - -a e ^ ^ ' 



-1 + 26e"" + 462e2(°^+2^'=")£;i(-26e''^) 



A, 



and 



(65) 



Pphoton = 3pphoton- (66) 

The resulting density of the dust reads, 

pdust = -2a2(6 + Xf^be"^ + ih^e^'''') + a2£,2g-2(az+26e") r _^ 4^paz _^ g^2g2(az+2fce''==)£;-(_2^gaz-)l _ 4^^^ (gy) 

where the constant D denotes the amplitude of the phantom, and D\ yields a effect of a cosmological constant, 
which can be understood as an AdS background. For a special configuration which can match Taub solution 6 = 
—36""^°, a = —1/3, zq = 1, it is not difficult to find positive definite densities for the dust and photon since we have 
two constants to adjust. For example, Di < 51.05 will ensure Pdust > in the interval z G [0, 1] when D = \. 



V. CONCLUSIONS AND DISCUSSIONS 

In this letter we explore the geodesic properties and the possible substances of the solution presented in [l| . 

We give the exact form of geodesies along z direction, both for the time-like and null cases. We find that the global 
effect of gravitation of the space put forward in [l[ not only the local effect, can be simulated by a field of a uniformly 
accelerating observer. 
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We discuss different cases of the source matter: the dust and photon, a scalar, a phantom, a phantom with dust 
and photon. We find that neither a single scalar nor a single phantom can play the role of source matter. We find 
that a phantom can sever as the source with dust and photon. 
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